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ABSTRACT 

The lattice data for the energy density of SU (2) gauge theory are calculated 
with non-perturbative derivatives of the coupling constants. These derivatives are 
obtained from two sources : i) a parametrization of the non-perturbative beta 
function in accord with the measured critical temperature and A/3— values and 
ii) a non-perturbative calculation of the presssure. We then perform a detailed 
finite size scaling analysis of the energy density near T c . It is shown that at the 
critical temperature the energy density is scaling as a function of VT 3 with the 
corresponding 3d Ising model critical exponents. The value of e(T c )/T^ in the 
continuum limit is estimated to be 0.256(23). In the high temperature regime the 
energy density is approaching its weak coupling limit from below, at T/T c pa 2 it 
has reached only about 70% of the limit. 
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1. Introduction 



The energy density for finite temperature SU(N) gauge systems is essential 
for the understanding and the investigation of the deconfinement transition and 
the plasma phase at high temperatures. In particular, a reliable estimate of the 
energy density e at the critical temperature T c is needed for the design of heavy 
ion collision experiments aiming at the production of the quark-gluon plasma. For 
theory e. g., the high temperature behaviour of e is of interest for the test of 
and the connection to finite temperature perturbation theory. The determination 
of the correct non-perturbative energy density is therefore an important goal of 
Monte Carlo simulations. 

The lattice calculations of e performed so far, though they are showing the 
same general behaviour for e on different size lattices, are not yet scaling, i. e. they 
do not lead to a unique function e(T). 

One of the sources of this disagreement has been addressed in [1] and is due 
to the use of weak coupling approximations to the derivatives with respect to 
the lattice anisotropy £ = a/a T ( a and a T are the lattice spacings in space and 
time directions ) of the space-like (<7 CT ) and time-like (g T ) coupling constants. The 
derivatives appear in the lattice expressions for the pressure P and the energy 
density 
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Here P„ yT are the space and time plaquettes, P CT]T = -^(Tr(l — U\U2U^U\)), 
on lattices of size N% x iV r . The plaquette Po is measured on an ^-lattice and 
serves to normalize the interaction measure A at T = 0. The derivatives have to 
be taken at £ = 1. Their sum is related to the f3- function through 

Once the (3— function and the pressure are known non-perturbatively, the en- 
ergy density can be calculated with the appropriate non-perturbative derivatives. 
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This has been done in [1] for SU(3) and JV^x4 lattices. In SU(2), where much 
more data on lattices with various N a and N T exist, this correction is still lack- 
ing. In the next two sections we shall therefore derive first a non-perturbative 
(3— function for SU{2) and then the non-perturbative derivatives separately. 

The other major effect on the energy density is coming from the finite size 
dependence on the spatial (N a ) and temporal (N T ) extents of the lattice. Here we 
have to distinguish two regions of temperature or coupling constants. In the neigh- 
bourhood of the deconfinement transition we expect finite size scaling governed by 
the critical exponents of the transition and the increase of the correlation length. 
In the weak coupling or high temperature regime we expect a different behaviour 
of the energy density, which is dominated by the weak coupling expansion of e and 
its N ff - and (stronger) iV T -dependence. We discuss the corresponding finite size 
equations in section 4. 

We shall compare the available SU (2) data for the energy density after cor- 
rection with non-perturbative coupling derivatives to the expectations from finite 
size scaling theory as described above. The results are presented and discussed in 
section 5. 

2. Beta function and scaling in SU(2) lattice gauge theory 

It is well known, that physical quantities such as the critical temperature or 
the string tension are violating asymptotic scaling in the range of coupling con- 
stants considered so far in lattice calculations. Nevertheless the scaling of dimen- 
sionless ratios of observables works remakably well [2] suggesting that deviations 
from asymptotic scaling may be described by a common non-perturbative (3- 
function. To parametrize the corresponding non-perturbative dependence of the 
lattice spacing a on g 2 we make the following ansatz 

aA L = R(g 2 ) ■ X(g 2 ) , (4) 

where 

^) = exp[-^ln(^)-^] (5) 

is resulting from the integration of the two-loop weak coupling expansion of the 
f3- function 
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-Pf = b g 3 + hg 5 + 



(6) 



Here 

UN , 34 iV 2 

are scheme independent. The inclusion of terms of order 0(g 7 ) in the expansion 
of the /3-function is then equivalent to a deviation of A from unity. 

The relation of the correction factor A to the sum of non-perturbative deriva- 
tives is found by rewriting Eq. (3) as 

and using 
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It is then convenient to write A for small g 2 in the following asymptotic form 

K9 2 ) = exp[^( Cl ^ 2 + c 2 g 4 + c 3 g 6 + ...)] , (10) 
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which implies A(0) = 1. 



The function A may be determined from any measured quantity of non-zero 
dimension. However only for small enough values of g 2 one may hope to find a 
unique function from different observables. 

In SU(2) lattice gauge theory we shall exploit the measurements of the critical 
coupling constant g 2 for iV T = 2,3,4,5,6,8,16 [2,3,4] to derive A. From the g 2 - 
values we find the critical temperature T c . On a lattice with iV r points in the 
temporal direction the temperature is 

Therefore we have 



M£) ^ ' (10) 
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Demanding independence of g 2 for T c /Al allows us then to determine the function 
X(g 2 ) in the presently studied coupling regime. In Fig. 1 we show the critical 
temperature 1/(N T R) as calculated from the two-loop formula, Eq. (5), i. e. by 
assuming A = 1, at the measured critical coupling for various N T . 

It has been noted [2,5,6] that the deviations from asymptotic scaling, which 
are obvious from Fig. 1 can, to a large extent, be accounted for by the introduction 
of an effective coupling constant. 

For our present analysis of bulk thermodynamics it is, however, more conve- 
nient to achieve a parametrization of the scaling violations in terms of the bare 
coupling constant. This has the advantage of giving us directly the (3— function, 
which we need anyhow for the evaluation of the energy density and the pressure. 

We thus proceed in the following way. First we use only those points, which 
belong to the weak coupling regime (4/g 2 > 2.23, i. e. N T > 4) to make a fit 
to 1/(N T R) with the asymptotic form, Eq. (10), of A. We have tried various 
possibilities for the exponent of A in Eq. (10). A rather good fit is obtained with 
the simple parametrization 

ci = c 2 = c n>3 = , (12) 

so that 

(13) 

The best result is given by C3 = 5.529(63) • 10~ 4 , which leads to 

T C /A L = 21.45(14) . (14) 

The last number is compatible with a corresponding extrapolation to a = in 
ref.[2]. 

The simple asymptotic fit Eq. (13) deviates in the neighbourhood of the 
crossover point from the measured two-loop critical temperature values, as can be 
seen in Fig. 1. In a second step we have therefore essentially fitted the measured 
points with a Spline interpolation, which continuously extends into the region, 
where Eq. (13) is valid. The solid line in Fig. 1 shows the final interpolation and 
the ratio of the measured l/(N T R(g 2 )) over the fitted X(g 2 ). The ratio should 
be a constant and equal to T C /A L . In Table 1 we give A and the derivative from 
Eq. (8), which up to a factor (-2) equals to the sum of derivatives. The latter is 
shown in Fig. 2. 
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Table 1 

The correction factor A and non-perturbative coupling derivatives. 

The interpolation was checked further by a comparison of the coupling constant 
shift A/3 to MCRG measurements with scale two blocking transformations [7,8]. 
In Fig. 3 the data and our A/3 are shown together with the two-loop prediction. 
At high (3- values the measurements from the small (8 4 ) lattice blocking are slightly 
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inconsistent with those on a large (32 4 ) lattice. Our result prefers the large lattice 
data and is, in contrast to the two-loop function, in full agreement with the rest of 
the data. In addition we have included in the figure data derived from the critical 
couplings g 2 of the deconfinement transition for iV r -values, which differ by a factor 
two. They are of course coinciding with our curve. 

3. The non-perturbative coupling derivatives 

The non-perturbative /9-function, which we have just derived, yields one of the 
two equations for the determination of the non-perturbative coupling derivatives. 
Another one is found by expressing, via Eqs. (1) and (2), the pressure in terms of 
the derivatives 
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In ref.[l] it was demonstrated how the pressure may be calculated as well non- 
perturbatively without derivatives, at least on large lattices. To this end one has 
to calculate the normalized free energy density from 
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where (3 = 4/g 2 , and to use the identity 



P=-f , 



(17) 



which is valid for homogeneous systems. It was checked for SU(2) [1] that a lattice 
of size 18 3 x 4 ( where we have many data points ) is large enough to justify this 
approximation. We have taken advantage of this fact and first calculated the 
pressure non-perturbatively on the 18 3 x 4 lattice by interpolating and integrating 
over [2P — (P a + P T )]. Subsequently we used Eqs. (8) and (15) to obtain 
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separately. It is clear, that below the critical point (4/g 2 < 2.30) , where both the 
pressure P and the plaquette difference P a — P T are small, the error on dg~ 2 /d^ 
becomes relatively large. We have therefore smoothly interpolated the values in 
that region. The resulting non-perturbative derivatives have been listed in Table 
1. We observe strong deviations from the asymptotic predictions [9] up to large 
values of the coupling (4/g 2 ~ 3.0). 

4. Finite size scaling behaviour 

A. The phase transition region 

On a finite lattice the correlation length is limited by the characteristic length 
scale of the system, i. e. normally the spatial extension N a . This limitation pre- 
vents the divergence of the correlation length at the critical point of a second 
order deconfinement transition. As a consequence the thermodynamic observables 
exhibit close to the transition finite size scaling properties, which are determined 
by the critical exponents and the length scale. According to the universality hy- 
pothesis of Svetitsky and Yaffe [10] the exponents of SU(2) should be the same as 
those of the three-dimensional Ising model. 

Usually finite size scaling analyses are carried out at fixed N T , but varying 
Na, see e. g. ref.[3]. The spatial extension determines then obviously the scale. In 
investigations, where results from lattices with varying N a and N T are compared, 
the dimensionless ratio 

^ = LT , (20) 

where L = N a a, has turned out to be the appropriate variable. This was observed 
already in studies of the heavy quark potential [11], and later on confirmed by the 
behaviour of the Binder cumulant [2]. 

We may check this idea further by looking at the interaction measure A, 
Eq. (2), for various N a and N T in the neighbourhood of the critical point, as a 
function of T/T c . Note, that here the non-perturbative /3-function enters in both 
the quantity A and the scale T/T c , but no single coupling derivative. In Fig. 4 
we see a clear coincidence of A-data belonging to the same ratio N a /N T . 

The finite size dependence of the energy density is related to the singular part 
f s of the free energy density. Since the critical exponent a of the specific heat is 
relatively small (a 0.11) [3], one expects that the regular parts in the specific 
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heat and in the energy density are dominating. Correspondingly, we make the 
following ansatz for the energy density close to T c 
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where we have already neglected irrelevant scaling fields. Here 



t=(T-T c )/T c 



(22) 



is the reduced temperature, Q e the scaling function of the energy density and v 
the critical exponent of the correlation length. 

B. The weak coupling region 

For small bare couplings g 2 , in the vicinity of the continuum limit, the energy 
density may be approximated by its weak coupling expansion. To obtain it up to 
order g 2 , we first derive e/T 4 from Eqs. (1) and (2) 
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Into this equation we insert the weak coupling expansions of the plaquettes for 
SU(N) lattice gauge theory [12] 

P = /^^ (2) +9\N 2 - !)P^) + ^ (2iV 2 -3)(iV 2 -l) p (46) + 0(/); (24) 

and the asymptotic values of the coupling derivatives [9] 
dg~ 2 dg~ 2 N 2 - 1 
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are related to those of the plaquettes in the following way 
a = 6(iV 2 - l)iV T 4 APi 2 r ) , 

a 1 = 6(iV 2 -1)N^NAP^ + 2 J^AP^ - 0.14671 1^^AP( 2 ) ^ 
+ AN [ 0. 000499 AP$ + 0.005306AP ( J ) ] 



where 

AP a/3 = P a -Pf3, (28) 

was used as an abbreviation and Po is the plaquette on the respective symmet- 
ric iV 4 lattice. The coefficients ao,a± approach with increasing N a and N T the 
corresponding continuum values 

«o=^(iV 2 -l) , 

15 (29) 
ai =--(N>-l)N. 

We have calculated both coefficients on a large number of N% x N T lattices. 
In Figs. 5 and 6 they are plotted for SU{2) at fixed ratios N a /N T as a function 
of N T . We observe, that ao has essentially the same ^-dependence for all ratios 
N a /N T and that the thermodynamic limit N a — > oo is almost reached already for 
N a /N T > 4. For a lattice calculation with an energy density as close as possible 
to the continuum limit one should therefore choose N T as large as possible at a 
moderate ratio N a /N T . 

Asymptotically, for N a — > oo and large N T the leading terms of ao are given 
for SU{N) by 



a = (N 2 - 1) 



7T 2 27T 4 1 / 1 \ 
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(30) 



The N a - and iV T -dependences of a\ are more complicated. For N a — > oo and large 
N T the second coefficient approaches the continuum value from below, in contrast 
to ao- For fixed N a /N T , both ao and a± cross their respective continuum limit 
lines, but at different iV r -values. 
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5. Numerical results for the energy density 



With our results for the non-perturbative derivatives and the /3-function we 
are now able to reevaluate the energy density with the plaquette values from refs. 
[3,4]. In Fig. 7 we show e/T 4 in the vicinity of the deconfinement transition as a 
function of T/T c . Like in the case of the interaction measure A we observe that 
at fixed temperature e/T depends only on the ratio N a /N T . From the ansatz, 
Eq. (21), the following scaling behaviour is then expected at t = 

£ = {*)A%) l °~ l)/ '- Q - lo) - (3l) 

We have made fits of this form to data at T c coming from seven different lattice 
sizes. Fixing the critical exponents to those of the 3d Ising model, i. e. 

(a-l)/u= -1.41 , (32) 

we find for the value of the energy density in the limit N a /N T — > oo ( the regular 
term at T c ) 

(j^j = 0.256(23) . (33) 

Our result is somewhat lower than the one found previously in ref. [13], where 
perturbative derivatives were used. A plot of the fit and the data are shown in 
Fig. 8. If we consider the exponent of (N a /N T ) in Eq. (31) as a free parameter 
we find as best value 

(a- l)/u = -1.58 . (34) 

This entails a slight change in the continuum limit value, Eq. (33) , which however 
was already taken into account in the error estimate. 

In the confined, low temperature phase an often used model for the energy 
density is that of a free gas of massive glueballs of various spins and parities. 
Assuming Boltzmann statistics the energy density of such a glueball gas is given 
by 
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In pure SU(2) lattice gauge theory the following glueball states have been 
identified : M Q+ /T c ~ 5.5, M 2 +/T c ~ 8.0, M Q -/T C ~ 9.5, and M 2 -/T c ~ 10.2 
[14] . At the critical temperature the energy density of an ideal gas of these glueballs 
is then found to be 

e gb (T c )/T* = 0.050 , (36) 

i. e. a factor five less than our estimate, Eq. (33). The presence of strong glueball 
interactions is therefore to be expected in the confinement phase. 

Comparing the critical energy density with that of an asymptotically free 
gluon gas 

e fl (T)/T 4 = 7r 2 /5, (37) 

we find a difference of nearly a factor eight. This means, that also in the decon- 
finement phase close to the transition strong interactions of the constituents of the 
gluon plasma are anticipated. 
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Table 2 

Results from Monte Carlo calculations on a 32 3 x 8 lattice. 

Let us now consider the high temperature phase. In Fig. 9a we show the approach 
of the energy density to the leading weak coupling term ao = e wc /T 4 - the lattice 
Stefan-Boltzmann limit - for two lattices. The 32 3 x 8 data are new. Here the 
lattice size was chosen such as to have a large iV T -value and a reasonable N a /N T - 
ratio at the same time. Details of these data are presented in Table 2. We see 
from Fig. 9a that the data even at T/T c w 2 are still about 30% below that of 
a gas of non-interacting SU(2) gluons and that the deviations from it at fixed 
temperature are only slightly iV T -dependent. Yet, we are still not in the weak 
coupling region, since the differences cannot be explained by the next-to-leading 
order weak coupling term. This does not come as a surprise, because in the T-range 
shown, the coupling derivatives are still non-perturbative. 
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In Fig. 9b we have plotted the difference between e wc /T 4 and the energy 
density data, showing again that the lattice Stefan-Boltzmann limit is clearly 
approached from below. A comparison of (e wc — e)/T 4 to finite temperature per- 
turbation theory with only one term proportional to the temperature dependent 
running coupling constant leads to a large value for g 2 (T) of the order of 4 ~ 5. 
It seems therefore that at presently accessible temperature values low order per- 
turbation theory is not applicable. 

6. Conclusions and discussion 

We have studied the finite size scaling behaviour of the energy density in 
SU (2) gauge theory at finite temperature. We find that the energy density at the 
critical point is only about 1/8 of the ideal gas value and also at T ~ 2T C it is still 
30% below this limit. The discrepancy between this result and earlier findings of 
a rapid approach to the ideal gas limit can be traced back to our non-perturbative 
definition of the pressure, Eq. (16). The additional use of a non-perturbative 
(3— function in the definition of (e — 3P)/T 4 is a minor modification to this. In 
fact the non-perturbative calculation of P/T 4 through Eq. (16) yields above T c a 
smaller value than Eq. (15) with perturbative coupling derivatives. This is reflected 
in a strong deviation of our results for dg~ 2 / <9£ from the perturbative values ( see 
Table 1). A similar result was obtained in SU(3), where these derivatives were 
directly calculated from Wilson loops on anisotropic lattices [15]. 

Some uncertainty on the absolute value of the pressure and consequently 
also the energy density is due to our approximation P/T 4 = for T ~ 0.8T C . 
This is suggested by the Monte Carlo data and seems also justified because the 
confinement phase consists only of heavy glueball states. We also note that unlike 
in the perturbative definition of thermodynamic quantities the non-perturbative 
approach does not automatically ensure that the ideal gas limit is reached for 
4/g 2 — > oo. However, it was checked in [1], that the order g 2 correction to Eq. (16) 
in the infinite temperature limit has the correct perturbative form. 

In the near future more Monte Carlo calculations in the high temperature 
region are needed to clarify the nature of the approach to the region where finite 
temperature perturbation theory is valid. 
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Figure Captions: 

1. The measured critical temperature 1/(N T R) for A = 1 ( circles) and T c / Al 

( squares) obtained after correction with the factor A ,vs. 4/g 2 for different 
N T . The solid line results from our interpolation of A, the dashed-dotted line 
from the asymptotic form A as . The dashed line is the fit result from Eq. (14). 

2. The sum of derivatives (dg~ 2 /d^ + dg~ 2 /d^) at £ = 1. The solid line is our 
non-perturbative result, the dashed-dotted line comes from A as and the dashed 
line from the two-loop approximation. 

3. The coupling constant shift A/3 from scale two blocking transformations vs. 
4/fir 2 , plotted at the larger of the two values. The data were measured by 
using two different lattices (diamonds, triangle) [7], blocking on one large 
lattice (crosses) [8] or obtained from T c (circles) [2]. The notation of the lines 
is the same as in Fig. 2. 

4. The interaction measure A = (e — 3T)/T 4 for various N a and N T in the 
neighbourhood of the critical point, as a function of T/T c . The necessary 
plaquette data were taken from refs.[3,4]. 

5. The first coefficient ao of the weak coupling expansion of e/T 4 versus N~ 2 
for the fixed ratios N a /N T = 2,4,8,16 (long dashes, dotted-dashed, dotted 
and solid lines). The continuum value is shown as a line of short dashes; the 
leading term in N~ 2 as a solid straight line. All numbers are for SU(2). 

6. The second coefficient a\ of the weak coupling expansion of e/T 4 versus N T 
for the fixed ratios N a /N T = 2, 4, 6, 8 (long dashes, dotted-dashed, solid and 
dotted lines). The continuum value is shown as a line of short dashes. All 
numbers are for SU(2). 

7. The energy density e/T 4 for various N a and N T in the neighbourhood of the 
critical point, as a function of T/T c . 

8. The energy density at T c as a function of (N a /N T )^ a ~ l ^ v with the critical 
exponents from the 3d Ising model. The straight line is a fit with the scaling 
form, Eq. (31). The notation for the data is like in Fig. 7, the diamond is 
from a 26 3 x 4 lattice. 

9. Comparison of the leading term of the weak coupling expansion ao = e wc /T 4 
to data for e/T 4 on two different lattices, versus T/T c . Fig. 9a shows the 
data directly, Fig. 9b the difference to the respective weak coupling limit. 
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